1/1 


AD-A133  10? 


COMPUTER  SIMULATION  OF  BUILDINGS  COOLED  BY  NATURAL 
VENTILATION(U)  NAVAL  CIVIL  ENGINEERING  LAB  PORT  HUENEME 
Cfi  S  ASHLEV  MAV  83  NCEL-TN-1666 


UNCLASSIFIED 


F/G  13/1 


NL 


urn;  mt  copy  nv~ n  13 3  /07 


N-1666 


t 


TN  NO: 


TITLE: 


COMPUTER  SIMULATION  OF  BUILDINGS 
COOLED  BY  NATURAL  VENTILATION 


AUTHOR:  S.  Ashley 
DATE:  May  1983 


SPONSOR:  Chief  of  Naval  Material 


PROGRAM  NO:  SO829-01-132A 


NAVAL  CIVIL  ENGINEERING  LABORATORY 
PORT  HUENEME,  CALIFORNIA  93043 


Approved  for  public  release;  distribution  unlimited. 


OTIC 

ELE'  r'- 

SEP  3  01983 


83  09  29 


030 


A 


METRIC  CONVERSION  FACTORS 


k1  in  >  2.54  (exactly).  For  other  exact  conversions  and  more  detailed  tables,  see  NBS 
Wise.  Pubi.  286,  Units  of  Weights  and  Measures,  Price  $2.25,  SD  Catalog  No.  Cl  3.10:286. 


Unclassified 

security  Classification  of  this  rage  'iwwi  n#»«  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  REPORT  NUM0FR  ]2  GOVT  ACCE£S»°N  NO. 

TN-1666  j  fiteSWff' '*1 

1  RECIPIENT’S  C  AT  AlOG  NUMBE  R 

4  Title  fmnd  Subtitle) 

COMPUTER  SIMULATION  OF  BUILDINGS  COOLED 

BY  NATURAL  VENTILATION 

S  TYPE  OF  REPORT  4  PERIOD  COVERED 

Not  final;  Oct  80  —  Sep  82 

6  PERFORMING  ORG  REPORT  NUMBER 

7  AUTHOR/"*. 

S.  Ashley 

6  CONTRACT  OR  GRANT  N  u  M  B  E  Rf  s  > 

9  PERFORMING  ORG*NI  Z  AT  ION  NAME  AND  ADDRESS 

NAVAL  CIVIL  ENGINEERING  LABORATORY 

Port  Hueneme,  California  93043 

to  program  element  project  task 

AREA  4  WORK  UNIT  NUMBERS 

63724N; 

SO829-01-1 32A 

It  CONTROLLING  OFFICE  NAME  AnO  ADDRESS 

Chief  of  Naval  Material 

Washington,  DC  20360 

12  REPORT  DATE 

May  1983 

»3  NUMBER  OF  PAGES 

53 

14  MONITORING  AGENCY  NAME  4  ADDRESS  fit  different  Irom  Controlling  Office) 

15  SECURITY  CLASS  t of  this  report) 

Unclassified 

IS*  DECLASSIFICATION  DOWNGRADING 

schedule 

16  DISTRIBUTION  STATEMENT  (of  thi*  Report) 

Approved  for  public  release;  distribution  unlimited. 

Aeces1 

11113 

*  T  .5 
a-  •  • 

l'1  '  d0  -  >  1 

17  DISTRIBUTION  STATEMENT  (of  the  mbatrset  entered  in  Block  20.  if  different  from  Report) 

J  . 

— L‘ .  a  1  .. 

18  Supplementary  notes  ; 

|  *  * 

19  KEY  WORDS  (Continue  on  reverie  aid •  it  nece**»ry  end  identity  bv  block  number )  ,  \ 

Natural  ventilation  cooling,  “stack  effect,”  wind  pressures  on  buildings,  effective  ■  \ 

temperature  limits 

20  ABSTRACT  (Continue  on  reverae  aide  If  nece*»*ry  end  identify  6v  block  number) 

'This  report  presents  general  information  covering  the  experimental  and  theoretical  basis 
used  to  develop  the  computer  simulation  of  natural  ventilation  cooling  of  buildings.  Computer 
input  and  output  and  an  application  example  are  included.  The  computer  program  predicts 
ventilation  rates,  mean  interior  air  speed  and  direction  of  air  flow,  interior  temperatures,  and 
inferior  and  exterior  effective  temperatures  for  a  given  set  of  conditions.  These  conditions 
are;  wind  velocity,  ambient  air  temperatures,  air  leakage  characteristics  of  openings,  the 

-  -  -  — continued - 

loll  F  T 

Ci  --  ’ 

* 


>  ;  •, 
a  t 

B  *  t 


3 


DO  1473  EDITION  OF  1  NOV  85  IS  OBSOLETE 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  flW»*n  Dst»  Entered) 


iMAttaJAlkJAAiJu  in  A 1  1 


jecumTr  classification  of  this  pagecwiab  dm*  i*i*t*d) 


20.  Continued 

-  mechanically  induced  air  flow  rates,  type  of  building  construction,  interior  heat  loads,  and 
the  number  of  people  in  the  building. 

The  application  of  this  computer  model  includes  the  study  of  the  natural  ventilation 
performance  of  a  building,  the  influence  of  natural  infiltration  on  mechanical  ventilation 
systems,  and  the  study  of  smoke  in  the  event  of  a  fire. , 


Library  Card 


Naval  Civil  Engineering  Laboratory 

COMPUTER  SIMULATION  OF  BUILDINGS  COOLED  BY  NATURAL 

VENTILATION,  by  S.  Ashley 

TN-1666  53  pp  iilus  May  1983  Unclassified 

1.  Natural  ventilation  cooling  2.  “Stack  effect”  I.  SO829-01-1 32A 

This  report  presents  general  information  covering  the  experimental  and  theoretical  basis  used 
to  develop  the  computer  simulation  of  natural  ventilation  cooling  of  building.  Computer  input 
and  output  and  an  application  example  are  included.  The  computer  program  predicts  ventilation 
rates,  mean  interior  air  speed  and  direction  of  air  flow,  interior  temperatures,  and  interior  and 
exterior  effective  temperatures  for  a  given  set  of  conditions.  These  conditions  are:  wind  velocity, 
ambient  air  temperatures,  air  leakage  characteristics  of  openings,  the  mechanically  induced  air  flow 
rates,  type  of  building  construction,  interior  heat  loads,  and  the  number  of  people  in  the  building. 

The  application  of  this  computer  model  includes  the  study  of  the  natural  ventilation  per¬ 
formance  of  a  building,  the  influence  of  natural  infiltration  on  mechanical  ventilation  systems,  and 
the  study  of  smoke  in  the  event  of  a  fire. 


Unclassified 


security  classification  of  this  ragE'hti «n  D Fntt'ctt) 


CONTENTS 

Page 


Chapter  1  INTRODUCTION  .  1 

Chapter  2  WEATHER  .  3 

2.1  WEATHER  DATA  ACQUISITION  .  3 

2.2  WIND  EFFECTS  ON  BUILDINGS .  3 

2.2.1  Nature  of  Wind .  3 

2.2.2  Wind  Pressures .  5 

2.2.3  Stack  Effect  .  6 

2.3  REFERENCES .  6 

Chapter  3  AIR  FLOW  CHARACTERISTICS  .  8 

3.1  EQUATIONS .  8 

3.2  REFERENCES .  9 

Chapter  4  EFFECTIVE  TEMPERATURE  .  12 

4.1  INTRODUCTION . 12 

4.2  OCCUPANT  METABOLIC  HEAT . 13 

4.3  BUILDING  HEAT  BALANCE . 15 

4.4  REFERENCES . 17 

Chapter  5  COMPUTER  SIMULATION  OF  NATURAL  VENTILATION 

COOLING  OF  BUILDINGS  .  20 

5.1  ASSUMPTIONS,  LIMITATIONS,  AND  BACKGROUND  PROGRAM  ....  20 

5.2  NETWORK  REPRESENTATION  OF  BUILDING  .  22 

5.3  REFERENCES . 23 

Chapter  6  COMPUTER  PROGRAM  .  25 

6.1  COMPUTER  PROGRAM  NOMENCLATURE  .  25 

6.2  ORDER  OF  INPUT  DATA . 27 

6.3  COMPUTER  INPUT  AND  OUTPUT  EXAMPLE  .  28 

6.4  CONCLUSION . 29 

APPENDIXES 

A  -  Weather  Data  Examples . A-l 

B  ■-  Pressure  Coefficients  for  Residential  Housing  and 

Commercial  Buildings  .  B-l 

C  -  Discharge  Coefficients  and  Air  Leakage  Data  for 

Building  Components  .  C-l 


v 


Chapter  1 
INTRODUCTION 


This  report  presents  information  covering  the  experimental  and 
theoretical  basis  used  to  develop  the  computer  simulation  of  natural 
ventilation  cooling  of  buildings.  Computer  input  and  output  and  an 
application  example  are  included. 

The  natural  ventilation  cooling  of  a  building  is  dependent  on 
(1)  weather,  (2)  shape  and  orientation  of  building,  and  (3)  shape  and 
position  of  doors  and  windows  (and  whether  open  or  closed).  Other 
influences  include  the  operation  of  mechanical  ventilation  systems  and 
the  transient  effects  caused  by  opening  and  closing  doors. 

Evaluation  of  rates  of  ventilation  and  interior  effective  tempera¬ 
tures  for  all  but  the  very  simplest  of  buildings  is  not  practicable 
without  the  use  of  computing  techniques. 

The  computer  program  predicts  ventilation  rates,  mean  interior  air 
speed  and  direction  of  air  flow,  interior  temperatures,  and  interior 
effective  temperatures  for  a  given  set  of  conditions.  These  conditions 
are  namely  meteorological  wind  speed  and  direction,  air  temperature,  air 
leakage  characteristics  of  openings,  the  mechanically  induced  air  flow 
rates,  type  of  building  construction,  and  the  number  of  people  in  the 
building. 

The  application  of  this  computer  model  (Figure  1.1)  includes  the 
study  of  the  natural  ventilation  performance  of  buildings,  the  influence 
of  natural  infiltration  on  mechanical  ventilation  systems,  and  the  study 
of  the  movement  of  smoke  in  the  event  of  a  fire. 
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Figure  I.].  Computer  program. 


Chapter  2 
WEATHER 


2.1  WEATHER  DATA  ACQUISITION 

Weather  data  information  on  surface  winds  and  psychrometric  summaries 
for  a  region  can  be  obtained  from: 

(1)  Local  weather  stations 

(2)  Naval  Weather  Service  Detachment,  Asheville,  N.C. 

Yearly  weather  data  are  derived  from  daily  observations  and  are 
presented  by  month  and  year  for  all  years  combined.  Examples  of  received 
data  are  shown  in  Appendix  A. 


2.2  WIND  EFFECTS  ON  BUILDINGS 

Wind  incidence  on  a  building  causes  positive  pressures  on  the  wind¬ 
ward  faces  and  negative  pressures  (or  suctions)  on  the  leeward  faces. 

As  a  result,  a  pressure  difference  is  generated  which  gives  rise  to  a 
basically  horizontal  movement  of  air  in  through  openings  or  cracks  in 
the  windward  side  of  the  building  and  out  through  similar  openings  in 
the  other  sides. 

2.2.1  Nature  of  Wind 

In  assessing  the  effects  of  wind  on  a  building,  it  is  important  to 
consider  the  characteristic  nature  of  the  wind.  The  wind  is  turbulent, 
and  its  mean  speed  varies  with  height.  Turbulent  characteristics  and 
the  vertical  profiles  of  wind  velocity  vary  with  the  stability  of  the 
atmosphere.  Local  topographical  features,  such  as  roughness  of  terrain 
over  which  the  wind  is  passing,  hills,  and  valleys,  can  also  affect  wind 
profiles.  Vertical  profiles  of  mean  wind  speed  for  boundary  layers  are 
approximated  by  taking  the  wind  speed  to  be  proportional  to  the  height 
raised  to  some  power  -  a  "power-law"  variation  (Ref  2.1).  The  simple 
expression  which  is  used  extensively  has  the  form 

U 

-jp  =  K  za  (2.1) 

m 
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where  =  wind  speed  from  local  meteorological  office 

=  wind  speed  as  specified  height 
K  =  terrain  factor 
a  =  terrain  exponent 
z  -  the  specified  height 


For  different  types  of  terrain,  Equation  2.1  describes  the  wind  speed 
variation  as  a  power-law  profile.  The  wind  speed  (U  )  is  measured  by 
the  local  meteorological  office,  and  is  measured  for  an  equivalent  height 
of  10  meters  (33  feet)  in  open  countryside. 

Equation  2.1  relates  the  wind  speed  at  any  other  height  (Uf)  and 
for  any  of  the  types  of  terrain  to  the  wind  speed,  U  .  Factors  for 
determining  mean  wind  speed  at  different  heights  and“1for  different  types 
of  terrain  are  given  in  Table  2.1. 

Table  2.1  Factors  for  Determining  Mean  Wind  Speed  at 
Different  Heights  and  for  Different  Types 
of  Terrain  From  Meteorological  Office  Wind 
Speed,  U  ,  Measured  at  10  Meters  in  Open 
Country 


Terrain 

K 

a 

Open  flat  country 

0.68 

0.17 

Country  with  scattered  windbreaks 

0.52 

0.20 

Urban 

0.  A0 

0.25 

City 

0.31 

0.33 

For  an  urban  site  and  a  building  height  of  20  feet,  the  wind  velocity  at 
the  roof  of  the  building  is  calculated  as  follows. 


U2Q  =  (0. AO) (20)°' 25  Um 

U20  =  °-85Um 


Hence,  the  roof  top  velocity  (or  local  velocity,  U  )  at  an  urban 
site  is  15%  less  than  the  measured  meteorological  wind  velocity  in  an 
open  site  for  the  same  gradient  velocity. 


A 


For  a  potential  building  site,  when  the  type  of  the  terrain  is 
established  (open  flat  country,  urban,  etc.)  the  appropriate  correspond 
ing  values  of  K,  and  desired  height,  z,  are  used  as  input  data  to  the 
program  and  the  expected  wind  speed  at  the  site  is  calculated  and  used 
to  calculate  wind  pressures  on  building  surfaces  and  expected  interior 
air  velocities. 


2.2.2  Wind  Pressures 

The  distribution  of  wind  pressure  around  a  building  depends  very 
closely  upon  the  local  variation  in  wind  velocity  which  the  building 
produces.  In  accordance  with  the  elementary  pressure-velocity  relation¬ 
ship  (Ref  2.2) 


Evidently,  as  the  velocity  (U)  increases  above  its  initial  value  (U  ), 
the  pressure  (P)  decreases  and  vice  versa  in  proportion  to  the  density 
(p)  and  the  square  of  the  velocity. 

For  convenience  the  pressure  distribution  around  the  building  is 
represented  by  a  dimensionless  pressure  coefficient  (C  ).  The  wind 
pressure  or  velocity  pressure  (Pw)  may,  therefore,  be  Represented  by, 


where  Pw  =  wind  pressure,  Pa  (pascals) 
p  =  air  density,  kg/m3 


(2.3) 


Ur  =  wind  speed  at  site  (usually  taken  at  roof  level),  m/sec 


Pressure  coefficients  applicable  to  simple  building  shapes  are  available 
(Ref  2.3).  For  residential  and  commercial  buildings  the  expected  pressure 
coefficients  are  given  for  different  wind  incidence  angles  in  Appendix  B 
and  References  2.1,  2.2,  2. A,  and  2.5. 

The  computer  program  will  calculate  the  wind  pressure  on  a  building 
surface  when  the  proper  pressure  coefficients  are  chosen  and  used  as 
input  data  to  the  program.  This  wind  pressure  in  turn  is  used  to  cal¬ 
culate  the  volume  flow  rates  in  each  building  room  (or  node) . 

To  determine  the  pressure  coefficients  applicable  to  the  surfaces 
of  complex  buildings,  it  is  necessary  to  use  wind  tunnel  testing 
techniques . 
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2.2.3  Stack  Effect 


The  difference  in  temperature  and  hence  in  density  between  air 
inside  a  building  and  air  outside  causes  a  movement  of  air  vertically 
through  the  building  via  such  openings  as  lift  shafts  and  stairwells. 
This  temperature-motivated  transfer  of  air  is  called  "stack  effects." 

The  pressure  differences  due  to  stack  effects  may  be  calculated, 
and  are  given  in  Reference  2.6  as 


(2.4) 


where 


P  =  pressure  difference,  Pa 


h  =  vertical  distance  between  inlet  and  outlet  openings,  m 


Tq  =  external  absolute  temperature,  °K 
T  =  internal  absolute  temperature,  °K 


If  the  temperature  difference  between  inside  and  outside  the  building  is 
20°C  and  the  story  height  is  3  meters,  then  the  pressure  difference  will 
be  2.6  Pascals  per  floor.  Hence,  for  a  10-story  building,  the  pressure 
that  will  be  added  to  the  wind  pressure  outside  is: 

9  x  2.6  Pa  to  all  ground  floor  wind  pressures 

8  x  2.6  Pa  to  first  floor  wind  pressures 

0  x  2.6  Pa  to  ninth  floor  wind  pressures 


For  most  natural  ventilation  cooling  of  building  applications,  the  "stack 
effect"  is  omitted.  The  "stack  effect"  is  included  in  the  computer  pro¬ 
gram,  and  it  is  an  option  to  the  user.  If  it  is  used,  then  the  pressures 
printed  at  the  end  of  the  program  will  be  essentially  incorrect;  however, 
the  air  volume  flow  rates  and  directions  of  flow  will  be  correct. 


2.3  REFERENCES 


2.1  A.  G.  Davenport.  1963  proceedings  of  the  Conference  on  Wind  Effects 
on  Buildings  and  Structures,  vol  1.  HMSO,  1965. 

2.2  S.  Goldstein.  Modern  developments  in  fluid  dynamics.  England,  Oxford 
University  Press,  1938. 


2.3  British  Standard  CP3,  Chapter  V,  part  2,  1972. 

2.4  University  of  Sydney.  Wind  generated  natural  ventilation  of  housing 
for  thermal  comfort  in  hot  climates,  by  R.  M.  Aynsley.  Sydney,  Australia. 
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2.5  National  Science  Foundation.  Grants  ENG72-04260-A01  and  ENG76-03035 
Wind  pressures  on  buildings,  by  R.  W.  Akins  and  J.  E.  Cermak.  Fall  1976. 

2.6  Department  of  Health  and  Social  Security.  BSRIA  Report:  Natural 
ventilation  in  hospital  buildings.  1976. 
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Chapter  3 

AIR  FLOW  CHARACTERISTICS 


3.1  EQUATIONS 

The  rate  of  air  flow  through  a  building  is  determined  by  the  air 
flow  from  high-pressure  to  low-pressure  regions  via  available  paths  and 
the  resistance  offered  by  the  flow  paths,  such  as  doors,  windows,  parti¬ 
tions,  etc. 

The  general  equation  relating  the  air  flow  (Q)  through  a  component 
with  a  pressure  differential  (AP)  across  it  (Ref  3.1)  is  expressed  by 
an  equation  of  the  form 

Q  =  K.  (AP)1/N  (3.1) 

where  Q  =  flow  rate 

Ki  =  leakage  coefficient  of  the  component 

AP  =  pressure  difference  between  two  rooms  (or  nodes) 

N  =  exponent  that  depends  on  the  nature  of  the  flow  through 
the  opening  and  ranges  between  1  and  2 


To  permit  use  of  this  generalized  equation  for  specific  situations, 
the  following  alternatives  may  be  considered. 

For  the  air  leakage  through  closed  doors  or  windows,  Equation  3.1 
is  often  expressed  as 


Q  =  C  l  (AP) 


(3.2) 


where  Q  =  flow  rate,  m3/sec 

C  =  crack  coefficient  per  unit  length  of  sash  crack, 
m3/sec/m  at  1  Pa 


£  =  crack  length,  m 


The  value  of  N  varies  between  1  and  2.  For  very  small  openings  (cracks), 
the  flow  is  laminar,  and  N  approaches  1.  As  the  size  of  the  opening  is 
increased,  the  flow  becomes  turbulent  and  N  approaches  2.  Values  of 
C  and  N  are  given  in  Appendix  C,  Table  C-2. 


For  open  windows  or  doors,  the  flow  rate  is  comparable  to  that 
through  a  thin-plate  orifice  with  flow  in  the  turbulent  region  and  a 
discharge  coefficient  of  0.65,  Equation  3.1  now  becomes 

Q  =  1.27  Cd  A.  (AP)1/2 


or 


Q  =  0.827  k±  (AP)1/2  (3.3) 

1/2 

where  0.827  is  the  leakage  coefficient  (m/s  -  (Pa)  ),  and  A.  is  the  area 
of  the  orifice  (m2).  For  different  values  of  discharge  coefficient,  C^, 
see  Appendix  C,  Table  C-2. 

For  air  flows  through  components  in  parallel,  as  shown  in  Figure  3.2, 
the  same  pressure  differential  acts  across  them  all,  and  the  total  air 
flow  rate  is  given  by 

Q  =  0.827  (AP)1/2  (Aj  +  A2  +  ...)  (3. A) 

where  A^ ,  k^,  etc.,  are  the  respective  orifice  areas  of  the  components. 

If  the  components  are  in  series,  the  effective  area  of  an  orifice 
in  Equation  3.3  is  given  by  (Ref  3.2) 


The  total  effective  area  of  the  orifice  is  therefore  always  less 
than  the  area  of  either  of  the  two  separate  components.  Figure  3.1 
(Ref  3.1)  shows  that  when  the  area  of  one  set  of  openings  is  more  than 
four  times  the  area  of  the  other  set,  then  only  the  area  of  the  smaller 
openings  needs  to  be  considered. 

Examples  of  how  to  combine  rooms/components  are  given  in  Figure  3.2 
(Ref  3.3).  In  Appendix  C,  Tables  C-l  (Ref  2.4)  and  C-2  (Ref  2.6)  give 
a  summary  of  air  leakage  data  for  building  components  taken  from  various 
sources. 


3.2  REFERENCES 

3.1  ASRAE.  Handbook  of  fundamentals,  chapter  25.  1967. 

3.2  B.R.E.D.  Principles  of  natural  ventilation.  Feb  1978. 

3.3  Department  of  Health  and  Social  Security.  BSRIA  Report:  Natural 
ventilation  in  hospital  buildings.  1976. 
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Option  A 


Standard  Method 


Conditions:  no  mechanical  ventila¬ 


tion  in  rooms 


from  node  1  to  node  3  A  =  Aj 
from  node  3  to  node  5  A  =  A2,  etc. 


from  node  1  to  node  3 


(Aj2  +  A22)1/2 


from  node  2  to  node  3 


(A22  +  A42)1/2 


Option  B 

Conditions:  Aj  *  Aj  and  A2  =  A4 

the  mechanical  ventilation  in  each 
room  is  equal;  wind  pressures  equal 


from  noe  1  to  node  2 
A  »  Aj  +  Aj 
from  node  2  to  node  3 
A  =  Aj  +  A A 


Option  C 

Conditions:  no  mechanical  ventilation 

in  rooms;  wind  pressures  equal 


from  node  1  to  node  2 


(Aj 2  +  A22)1/2  (Aj2  v  A42>' 


Figure  3.2.  Examples  of  how  to  combine  rooms  or  components. 


Chapter  A 

EFFECTIVE  TEMPERATURE 


A. 1  INTRODUCTION 

A  physiological  measurement  combining  the  effects  of  temperature, 
humidity,  and  air  movement  is  called  the  effective  temperature  (ET) 
scale.  New  comfort  limits  by  ASHRAE  Standard  55-7A  are  for  ET*  between 
72°F  and  78°F  with  maximum  relative  humidity  of  80%  for  the  upper  limit 
(Ref  A. 1) . 

The  effective  temperature  that  results  inside  a  building  depends  on 
many  factors,  including: 

1.  Ventilation  air  flow  rate  (Q)  through  the  building 

2.  Dry-bulb  temperature  (T^g)  and  wet-bulb  temperature  (T^g)  of 
ambient  air 

3.  Metabolic  heat  generation  rate  of  occupants  (Qqcc) 

4.  Heat  transferred  through  the  envelope  of  the  structure 

5.  Interior  heat  loads  (Q  ),  such  as  direct  sunlight,  electric 
light  bulbs,  etc. 

6.  Air  flow  pattern  inside  the  building,  etc. 

Heat  transferred  through  the  building  envelope  is  assumed  steady 
state.  Air  within  the  building  is  assumed  to  be  uniformly  distributed 
such  that  no  air  stratification  or  temperature  variation  occurs.  The 
building  effective  temperature  is  then  assumed  to  be  a  function  of  the 
following: 

EI  *  «V  %CC-  «•  TDB’  W  (41> 


Using  the  building  configuration,  the  hourly  data,  and  the  number  of 
occupants,  an  example  of  the  average  daily  effective  temperature  inside 
each  room  of  the  building  due  to  varying  ambient  weather  is  plotted 
against  time  as  shown  in  Figure  A.l. 
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4.2  OCCUPANT  METABOLIC  HEAT 

The  average  metabolic  rates  for  sedentary  activity  levels  as  func¬ 
tions  of  age  and  sex  for  the  most  recent  census  of  the  U.S  population 
are  given  in  the  1978  ASHRAE  Applications  Handbook;  Table  4.1  and  subse¬ 
quent  equations  were  developed  by  GATEX  (Ref  4.2)  using  these  data. 

Table  4.1.  Sedentary  Metabolic  Rates  as  a 
Function  of  Age  and  Sex 


Age 

Percent  of 

Population 

_ 

Metabolic 

Heat  (Btu/hr) 

Male 

Female 

Male 

Female 

75  and  over 

1.5 

2.2 

310 

70  to  74 

1.1 

1.5 

320 

65  to  69 

1.5 

1.8 

350 

60  to  64 

2 

2.3 

360 

55  to  59 

2.3 

2.7 

370 

240 

50  to  54 

2.6 

2.8 

375 

245 

45  to  49 

2.9 

3.1 

380 

250 

40  to  44 

2.9 

3 

395 

250 

35  to  39 

2.7 

2.9 

400 

250 

30  to  34 

2.8 

2.9 

400 

250 

25  to  29 

3.3 

3.3 

400 

250 

20  to  24 

3.9 

4.2 

405 

250 

15  to  19 

4.7 

4.6 

400 

260 

10  to  14 

5.3 

5 

275 

250 

5  to  9 

5 

4.8 

200 

200 

Under  5 

4.3 

4.1 

150 

150 

Total 

48.8 

51.2 

-- 

-- 

a"Social  Indicators,"  1973,  U.S.  Dept,  of  Commerce 


The  polynomial  expressions  for  metabolic  heat  generation  (Ref  4.3) 
are  as  follows: 


Qocc, total  =  a  (ET>  +  b  <ET>  +  c 


Qocc, sensible  =  A  ^DBB^  +  B  ^DBB^  +  C 


(4.2) 


where  T^gg  =  the  dry  bulb  temperature  of  the  building  air 
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1.  a  =  0;  b  =  1.482;  c  =  514,  for  50°  £  ET  £  87°F 

2.  a  =  -1.508;  b  =  259.7;  c  =  -10,795.2,  for  87°  <  ET  i  102°F 

3.  a  =  0;  b  =  0;  c  =  0,  for  ET  >  102°F 

4.  A  =  -0.06875;  B  =  1.625;  C  =  523,  for  Tpgg  £  50°F 

and 

Q  ,  =  0  +  0 

occ, total  occ, sensible  ^occ, latent 


Equations  4.2  and  4.3  have  been  applied  and  proved  for  only  healthy 
male  subjects  in  the  age  group  of  18  to  24  years.  For  this  study  an 
"average  person,"  whose  metabolic  heat  rate  is  the  weight  average  that 
has  been  created  using  the  following  equation: 


Q 


average  person 


(%  of  males)  x  (Btu/hr  male)  + 
(%  of  females)  x  (Btu/hr  female) 
100 


=  281.32  Btu/hr 


(4.4) 


This  is  70%  of  the  value  given  in  the  table  for  males  between  20  and 
24  years.  This  weight  average  is  used  in  the  program  to  give  the  percent 
of  "average"  people  who  will  be  comfortable  within  prescribed  tempera¬ 
ture  limits.  Using  the  above  equations  and  data  from  Ref  4.4  (pages  6 
to  8,  Tables  2  and  3)  the  following  equations  were  generated  for  calcu¬ 
lating  occupant  metabolic  heat: 

«„cc, total  *  0.413|a(ET)2  *  MET)  +  c]  IN 
+  1.091  T  u.0,6 

i=l  1  J 

(4.5) 

^oce, sensible  *  0  <'13'fl(TDBB»2  *  B(W  *  C'  '» 

+  1.091  T  u.0-6 

i=l  1  J 

where  N  =  number  of  occupants 

=  air  velocity  felt  by  occupant  i 


14 


The  effective  temperature  (ET)  of  the  building  air  can  be  calculated 
as 


ET  = 


107.5(TjjbB) 


45,2^twbb^ 


TWBB  +  62,3 


(4.6) 


where  ^jjgg  =  dry-bulb  temperature  of  building  air 


T^bb  =  wet-bulb  temperature  of  building  air 


4.3  BUILDING  HEAT  BALANCE 

Sensible  and  latent  heat  loads  and  ventilating  air  are  introduced 
into  the  building,  air  is  exhausted  out  of  the  building  and  energy 
is  lost  (see  Figure  4.2). 

The  following  assumptions  are  made  (Ref  4.2): 

1.  The  air  within  the  building  is  completely  and  instantaneously 
mixed . 

2.  The  film  heat  transfer  coefficients  are  constant  for  any  one 
boundary  surface. 

3.  The  radiative  energy  transfer  within  the  building  can  be 
neglected. 

4.  The  condition  of  the  air  exhausted  from  the  building  is  the 
condition  of  the  building  atmosphere. 

5.  The  thermal  and  physical  properties  of  the  structural  materials 
and  of  air  are  not  temperature-dependent. 

6.  The  incident  solar  radiation  is  absorbed  on  the  outer  surface 
of  the  boundaries  and  appears  as  conducted  energy  or  is 
transmitted  into  the  building  and  is  considered  an 
instantaneous  load  along  with  other  internal  loads. 

7.  The  thermal  loads  and  the  psychrometric  states  of  the  inlet 
and  exhaust  air  are  constant  over  short  time  intervals. 

A  heat  balance  of  the  building  volume  yields  the  following: 
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Sensible  Heat  Gain  =  Qocc>sensible  +  x  Fw  +  KR  x  FR) 

(tdb  ’  tdbb^ 

Latent  Heat  Gain  =  Q  ,  . 

occ, latent 

Sensible  Heat  Loss  =  60  x  x  p  x  Q  x  ^DB  ~  ^DBB^  (A .7) 

Latent  Heat  Loss  =  0.625  xlO^xpxQx  (WR  -  W) 

density  of  leaving  air 
flow  rate 

heat  transfer  coefficient  for  building  wall  and 
building  roof 

total  external  wall  or  roof  surface  area 
humidity  ratio  of  building  air 
humidity  ratio  of  outside  air 
specific  heat  of  air 


where  p  = 

Q  = 

kw&kr  = 

FW  &  FR  = 

WB  = 
W  = 


Since  we  assume  that  the  building  is  in  thermal  equilibrium  at  the 
end  of  each  time  interval, 

Heat  Gain  =  Heat  Loss  (4.8) 
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Once  the  flow  rate,  Q,  is  established,  the  resulting  interior 
effective  temperature  can  be  calculated  at  each  room  and  the  percent  of 
people  who  will  feel  comfortable  is  given  along  with  the  percent  of 
relative  humidity. 


4.4  REFERENCES 

4.1  ASHRAE.  Handbook  of  Fundamentals.  1977,  pp  8.1-8.36. 

4.2  GATEX.  CARD  Report  Al-11  (1713):  Adequacy  of  wind  ventilation  in 
upgraded  shelters.  Jun  1980. 

4.3  F.  C.  Houghten.  "Heat  and  moisture  from  the  human  body  and  their 
relation  to  air  conditioning  problems,"  ASHVE,  Transaction,  vol  35,  1929. 

4.4  ASHRAE  Handbook  of  Fundamentals,  Chapter  8,  Physiological  principles, 
comfort  and  health.  1977. 


Time  (hours  and  months) 


Figure  4.1.  Variation  of  building  (one  room)  daily  average  effective 
temperature  with  time. 


Chapter  5 


COMPUTER  SIMULATION  OF  NATURAL  VENTILATION 
COOLING  OF  BUILDINGS 


5.1  ASSUMPTIONS,  LIMITATIONS ,  AND  BACKGROUND  PROGRAM 

A  brief  summary  of  assumptions,  limitations,  and  formulation  of 
the  "CRKFLO"  computer  program  (Ref  5.1)  is  presented  in  this  chapter. 

1.  Net  air  supplied  (or  extracted)  by  mechanical  ventilation  system 
is  assumed  to  be  constant  and  independent  of  building  pressures. 

2.  Each  room  (or  node)  of  the  building  is  considered  to  be  at  a  uni¬ 
form  pressure  throughout  the  room. 

3.  If  very  large  openings  are  considered,  difficulty  of  convergence 
may  be  experienced.  This  restriction  can  be  overcome  by  redesigning  the 
node  network.  This  problem  occurs  if  pressure  of  less  than  10  Pascal 
is  encountered  across  the  component,  so  the  resultant  flow  will  probably 
be  insignificant. 

4.  The  average  interior  air  speed  is  calculated  by  using  air  flow 
rate  divided  by  cross-sectional  area  of  room  (node). 

5.  The  average  dry-bulb  and  wet-bulb  temperature  of  each  node  is  a 
function  of  incoming  air  conditions  at  the  room  and  represents  conditions 
of  air  leaving  the  room  (node) . 

6.  There  is  no  heat  transfer  through  partitions. 

A  network  of  air  flow  paths  comprises  the  nodes  and  their  associated 
interconnections,  which  in  many  ways  is  analogous  to  flow  through  electric 
circuits  (Ref  5.1).  The  solutions  of  these  are  in  general  complicated  and 
require  trial  solutions  in  which  the  elementary  circuits  are  balanced  in 
turn  until  all  conditions  for  the  flow  are  satisfied.  Three  basic  con¬ 
ditions  must  be  satisfied  in  any  network: 

1.  The  algebraic  sum  of  the  pressure  drops  around  each  circuit 
must  be  zero. 

2.  Flow  into  each  node  must  equal  the  flow  out  of  the  node. 

3.  The  equation  relating  the  pressure  difference  and  flow  for 
each  interconnection  must  be  satisfied. 
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The  solution  to  the  equation, 


Q  =  K  (AP) 


is  obtained  by  an  iterative  process  (Ref  5.2),  in  which  during  each  sue 
cessive  cycle  corrections  are  made  to  the  values  of  the  pressure  heads 
at  the  nodes  (or  rooms)  for  which  pressures  were  not  specified.  The 
source  of  the  corrective  term  is  as  follows: 

For  any  interconnection 

AP  =  APQ  +  CR  ( 

where  AP  =  correct  pressure  difference 
APq  =  assumed  pressure  difference 
CR  =  correction 

Then,  for  each  interconnection 


Q  =  K  (APq  +  CR) 


Q  =  K  (APq1/N  +~AP0(1/N  •  CR  + 


By  expanding  the  Taylor  series  of  Equation  5.3  to  include  the  second 
order  differential,  the  equation  becomes 

q  =  k  •  ♦£aP0(,/N  -  ■>  •  cr 

.  O/N  -  2)  .  „„2 


For  all  the  interconnections  to  a  room  (node), 

Q  =  IK  APq1/N  =  IK  APq1/N  +  CR  -fiAP0(1/N  ' 

+  CR2  •  #  7  V  ~  AP (1/N  -  2)  =  0  ( 


CR  has  been  taken  out  of  the  summation  since  it  is  the  same  for  all  the 
interconnections  around  the  room.  The  quadratic  solution  for  CR  becomes: 


AP0(1/N-1)  ±|/-jp  AP0(1/N-1)2  -  4  •  (IK  APq1/N)  •  (|^(*-1)  AP(1/N"2) 

o  .  O/M-l)  .  ApO/N-2) 

*  N  2 

. (5.6) 


The  use  of  this  correction  factor  yields  a  fast  convergence  to  the 
solution. 


5.2  NETWORK  REPRESENTATION  OF  BUILDING 

A  plan  view  of  the  building  under  consideration  is  obtained,  and 
each  room  or  compartment  represents  a  node.  Outside  wind  pressure  areas 
are  also  represented  by  nodes  (Ref  5.1).  Each  set  of  nodes  is  intercon¬ 
nected  by  air  flow  paths  (see  Figure  5.1).  It  is  assumed  that  substan¬ 
tially  uniform  pressure  conditions  prevail  at  each  node.  The  intercon¬ 
nections  correspond  to  impedance  to  the  air  flow  such  as  windows  and 
doors.  A  knowledge  of  the  air  flow  characteristics  is  required  for  each 
door  and  window  contained  in  the  network. 

The  number  of  exterior  nodes  is  designated  as  MH,  and  MT-MH  no.’.s 
represent  rooms  and  corridors  inside  the  building.  Pressures  for  the 
external  nodes  must  be  specified  and  these  pressure"  will  represent  a 
combination  of  wind  pressures  and  buoyancy  pressuvrr  « stack  »Hect). 

The  nodes  are  numbered  from  1  in  ascending  &r4er,  starting  with  the 
nodes  for  which  a  known  pressure  will  be  specified  (see  Figure  5.1); 
there  will  be  MH  of  these.  The  numbering  for  the  internal  nodes  (rooms 
and  corridors)  continues  in  sequence,  starting  at  MH+1  through  to  MT, 
where  MT  is  the  last  node. 

R  denotes  the  number  of  each  node  connected  to  P  taken  in  numeri¬ 
cally  ascending  order,  and  is  only  defined  if  R>P;  otherwise,  it  is  zero 
(see  Figure  5.1).  K  and  N  are  the  characteristic  constants  for  each 
particular  interconnection.  Zero  terminates  the  sequence  of  intercon¬ 
nections  to  each  particular  node  P  and  must  appear  even  if  no  inter¬ 
connections  are  listed. 

To  convey  the  network  shape  to  the  computer,  the  data  must  be  input 
in  a  strict  order.  The  following  is  an  example  of  the  generation  of  a 
network  of  Figure  5.1. 

Listing  Notes 

P  =  1 

R  =  4  R  >  P 

K14  N14 

R  =  0  R  =  0  terminates  the  sequence. 
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For  input  area,  A,  the  value  of  K  is 
computed  by  the  program. 


MH  =  1,  2,  and  3  are  the  exterior 
nodes  with  known  pressures. 


From  node  4  to  node  1,  R  is  less  than 
P  (1  <  4),  thus  R  =  0.  This  is  to 
avoid  listing  interconnections  twice. 


MH+1  =4,  5,  and  6  are  the  interior 
nodes  with  the  unknown  pressures. 


Pressures  are  to  be  specified  at  nodes  1,  2,  and  3. 


When  mechanical  ventilation  is  present  in  a  room  the  volume  flow 
rate  in  m3/s  is  specified  and  used  as  input  data  to  the  computer  program. 
Additional  input  data  consists  of  room  volume  and  length.  The  program 
calculates  pressures  and  pressure  differences  at  and  between  nodes, 
volume  flow  rate,  and  direction  air  changes  per  hour  at  each  node, 
and  interior  air  speed.  The  study  of  direction  of  smoke  in  the  event 
of  a  fire  in  the  room  can  be  predicted  by  the  direction  of  the  flow  in 
or  out  of  the  room. 


5 . 3  REFERENCES 

5.1  Department  of  Health  and  Social  Security.  BSRIA  Report:  Natural 
ventilation  in  hospital  buildings.  1976. 

5.2  G.  T.  Tamura.  "Computer  analysis  of  smoke  movement  in  tall  buildings 
ASHRAE  Transactions,  1969,  75,  II,  81. 
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Chapter  6 
COMPUTER  PROGRAM 

This  program  is  stored  at  CDC  KWA,  Sunnyvale,  Calif.,  on  PUBLIC 
READ  TAPE:  VSN  =  PN2339,  and  it  is  written  in  FORTRAN  IV  computer 
language.  The  program  consists  of  two  main  routines;  Routine  Flow  cal¬ 
culates  air  flow  through  the  building,  and  Routine  Air  estimates  tempera¬ 
tures  based  on  results  of  Flow  routine. 

6.1  COMPUTER  PROGRAM  NOMENCLATURE 

The  following  explains  each  data  element  in  order. 

AKW  Total  heat  transfer  coefficient  of  exterior  wall,  Btu/hr-ft2 

AKR  Total  heat  transfer  coefficient  of  roof,  Btu/hr-ft2 

ALF  Exponent  for  calculating  local  velocity 

AR  Area  of  inlets,  ft2  (if  "stack  effect"  used,  divide  calculated 

area  by  0.827) 

CP  Pressure  coefficient 

CQ  Flow  equation  exponent 

DBT  Dry-bulb  outside  temperature,  °F 

EL  Length  of  room,  ft 

FN  Number  of  specific  floor 

FR  Roof  area,  ft2 

FW  Exterior  wall  area,  ft2 

HB  Height  of  building,  ft 

HS  Initial  estimate  of  pressure  at  each  internal  node,  Pascal 

ITITL  Name  of  case 
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ITLIM 


Maximum  number  of  matrix  inversions  allowed 


KM  Terrain  constant 

M  Flag,  equal  to  zero  means  no  "stack  effect";  equal  to  one 

means  "stack  effect" 

MH  Number  of  nodes  at  which  pressure  is  specified 

MP  Maximum  number  of  interconnections  at  any  node 

MT  Number  of  nodes 

NAR  Number  of  interconnections  for  flow  equation  constants,  <  21 

NCS  Maximum  number  of  iterations  allowed 

NF  Number  of  floors,  including  ground  floor 

NUM  Number  of  occupants 

P  Node  (room) 

PB  Barometric  pressure  x  100,  in.  Hg 

Q  Mechanical  ventilation  at  each  internal  node,  m3/sec 

QL  Electric  heat  load,  kW 

QSOL  Algebraic  sum  of  other  heat  loads,  Btu/hr 

R  Node  (to) 

TI  Absolute  interior  temperature,  °F,  (for  M  =  1,  "stack  effect") 

TO  Absolute  exterior  temperature,  °F,  (for  M  =  1,  "stack  effect") 

ULT  Convergence  criterion  for  sum  of  flow  errors 

UM  Meteorological  wind  speed,  knots  or  mph 

VOL  Room  volume,  ft3 

WBT  Outside  wet-bulb  temperature,  °F 

WD  Wind  directions  in  degrees 

WFLAG  =  0  meteorological  wind  speed  in  mph 

=  1.  meteorological  wind  speed  in  knots 
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6.2  ORDER  OF  INPUT  DATA 
DECK  SET-UP 

NAME, PRIORITY, TIME, STK  WA. 

USER, USER  NAME, PASSWORD, KVA. 

CHARGE , CHARGE# , PROJECT# . 

COPYCR, INPUT, TAPES. 

COPYCR, INPUT, TAPE7. 

REQUEST , TAPE , VSN=PN2  339 , D=PE , UN=L6023GS . 

(JCL)  COPYBF , TAPE , PROGR . 

REWIND, PROGR. 

UNLOAD, TAPE. 

FTN ,R=3 ,L=0 , I=PROGR . 

ATTACH , CELLIB/UN=L5 1 07GS . 

LDSET , PRESET=ZERO , L I B=CELLIB . 

LGO,PL= 100000. 

ROUTE, PLOT ,UN=L0601GS ,ST=WCZ ,DC=PU . 

/*EOR 

(TAPE5)  WEATHER  DATA  (HOUR, DRY  BLB.WET  BLB, BAROM. PRESS. , MONTH) 

ON  LAST  RECORD  HOUR=99 
FORMAT  (15, 3F  10.1.1XA3). 

/*EOR 

(TAPE7 )  WIND  SPEED  (ONE  VALUE/RECORD) 

DATA  ON  LAST  RECORD  UM=999 

/*EOR 

TITLE  (§80  CHARACTERS)  IF  EMPTY=END  OF  JOB 
NF.FN 

HB,ALF,KM,M 

NAR 

(AR(I),I=1,NAR) 

MT,MP,HS 

P 

R  REPEAT  AS  OFTEN  AS  NECESSARY 

CQ  IF  R  *  0 

MH.ULT 

NCS,ITLIM 


Q(I) 

IF  M  *  0 
TI,TO 
I=MH+1 ,MT 

EL(I) 

I=MH+1 ,MT 

QL(I) 

I=MH+1,MT 

QSOL(I) 

I=MH+1,MT 

AKW(I) 

I=MH+1 ,MT 

AKR(I) 

I=MH+1,MT 

FW(I) 

I=MH+1,MT 

FR(I) 

I=MH+1 ,MT 

VOL(I) 

I=MH+1,MT 

NUM(I) 

I=MH+1 ,MT 

WFLAG 
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WD, (CP (I) ,1=1 ,MH)  REPEAT  THIS  RECORD  FOR  WIND  DIRECTION. 

LAST  RECORD  WD=999 
REPEAT  FROM  TITLE  RECORD  ON 

/*EOR 

END  OF  FILE=SPECIAL  CARD  IF  CARD-INPUT 
=/*EOR  IF  TERMINAL  TYPE  IN. 


6.3  COMPUTER  INPUT  AND  OUTPUT  EXAMPLE 

Using  the  network  representation  of  a  single  story  on-grs.de 
(Figure  5.1)  building  at  a  zero  degree  incidence  to  long  walls ,  (Ref  2.4) 
Appendix  B,  Table  B-2,  the  input  and  output  data  are  as  follows: 

EXAMPLE  VENTILATION,  ONE  STORY  BLDG,  GAMA  0.5 


11  CARD  LI S TIN 6  FRI,  SEP  24  1982  AT  13:15:41 

2  1  REMOTE  C  DC 
31  SOPMIAtP4»T?0O*S  TRIM. 

4)  USER  ,L60  23GS  ,NCE  LPAS  ,KMA. 

5)  CHARGE  ,M4821DS, 31463  0897  . 

S)  COPT CRfI NPUT « TAP E5. 

71  COPY  CR,I  NPUT  »  TAP  E  7. 

81  REQUEST, TAP:  , VSN =PN? 3 39, PC=W ,0=PE. 

91  COPY  BF,T  APE,  SOPH  I  A. 

10)  RE  HI  NO  ,S  OPHI  A  • 

Ill  UNLO  AD  ,T  APE. 
l?>  FTN,  R-  3,  I-S3PHIA  . 

131  ATTA  CH,C  £LLi  fi/UN=L5l  0  7GS. 

14  1  LOSE  1,  PRESET  =  7E3  0,U  R=CEL  LI  8. 

151  LGO, PL  =1  0003  0  . 

161  ROUT  c ,PL  OT ,JN  =L0  E01GS*ST  =  WC7 »OC  =  PU  • 


17) 

EXIT  . 

18) 

FXIT  . 

19) 

FTN.  P=  5,  I  ~  $3 

PMIA  ,  L 

=1  • 

20  ) 

/•  E3  R 

21) 

1  7'. 7 

70. > 

29.  °2 

JLL 

22  ) 

4  7P.2 

70.) 

27.  c? 

JLL 

23) 

7  7  £  .  0 

70.4 

29.  op 

JLL 

24  ) 

10  7  *:  •  0 

71 

29.  4? 

JLL 

25) 
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In  the  computer  program  output  the  direction  of  flow  is  indicated 
by  the  sign.  A  negative  sign  in  FLOW  indicates  that  air  flow  INTO  the 
room  (or  node).  Dry-bulb,  wet-bulb,  and  effective  temperature  for  air 
LEAVING  the  room  is  also  given. 

6.4  CONCLUSION 

This  computer  simulation  of  natural  ventilation  cooling  of  buildings 
can  be  used  as  a  tool  for  a  preliminary  study  of  a  building  design, 
to  determine  size  of  openings,  required  air  flow,  air  changes,  interior 
air  velocity,  and  comfort  ranges.  In  FY83  results  of  field  testing  and 
wind  tunnel  testing  will  be  compared  to  results  obtained  from  this 
compute r  program. 
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Table  A-l.  Means  and  Standard  Deviations 
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Table  A-2.  Means  and  Standard  Deviations 
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Appendix  B 

PRESSURE  COEFFICIENTS  FOR  RESIDENTIAL  HOUSING 
AND  COMMERCIAL  BUILDINGS 


Table  B-l.  Wind  at  45°  Incidence  to  a  Long  Wall 

(When  data  from  Table  B-l  and  Table  B-2  are  used,  set  the  pressure 
difference  coefficient,  ACp,  equal  to  the  windward  pressure  coefficient, 
Cpw,  and  the  leeward  pressure  coefficient,  Cp^  equal  to  zero.  So 

ACp  =  Cpw 

and 

CpL  =  0) 


Model  Features 

Average  Pressure  Differential 
Coefficient  Between  Mid-Height  of 
Long  Walls  for  Wind  at  45  degrees 
Incidence  to  a  Long  Wall  (Refer¬ 
ence  Pressure  -  Dynamic  10  Meters 
Above  Ground  in  Wooded  Suburban 
Terrain) 

Single  story  elevated  above  ground 
with  extended  eaves  and  end  walls 
(Type  4)a 

0.75 

Single  story  on-grade  with  extended 
eaves  and  end  walls  (Type  2) 

0.58 

Single  story  elevated  above  ground 
(Type  3)a 

0.43 

Single  story  on-grade  with  flat 
roof  (Type  6) 

0.35 

Single  story  on-grade  with  10 
degree  pitch  roof  (Type  1) 

0.34 

Two  story  (Type  5) 

0.31 

a 

Note  that  the  three  highest  ranking  model  types,  4,  2  and  3,  at  45 
degrees  incidence  were  also  highest  ranking  averaged  over  wind  inci¬ 
dences  of  0  degrees,  30  degrees  and  45  degrees. 
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2. 1  Pressure  Coefficients  for  Buildings  with  Three  Floors  or  Less. 


Test  data  for  side  ratios  of  y  =  1.0,  0.5  and  0.25  of  model  flat- 

roofed  rectangular  buildings  have  been  accumulated  and  given  in  Graphs  A, 

B,  and  C.  When  side  ratio  is  close  to  one  of  these  ratios,  select  the 
data  for  that  ratio. 

1.  Wind  direction,  angle  a.  Select  wind  incidence  on  long  side, 
side  3,  windward,  and  side  1  leeward  from  Figure  B-2b. 

2.  The  pressure  coefficient,  Cp^  windward  (Cp^)  is  read  from 
Graphs  A,  B,  and  C  depending  on  the  value  of  side  ratio.  The  pressure 
coefficient  Cp^  leeward  (Cp£)  is  taken  to  be  zero,  Cpj  =  0. 

3.  Wind  speed  is  calculated  at  roof  height. 

2.2  Pressure  Coefficients  for  Tall  Buildings. 

Test  data  for  side  ratios  of  1.0,  0.5  and  0.25  of  model  flat- 

roofed  rectangular  buildings  have  been  accumulated  and  are  given  in 

Graphs  D,  E,  and  F.  When  side  ratio  is  close  to  one  of  these  ratios, 
select  the  data  for  that  ratio. 

1.  Wind  direction,  angle  a.  Select  wind  incidence  on  long  side, 
side  3,  windward,  and  side  1  leeward  from  Figure  B-2b. 

2.  The  pressure  coefficient,  Cp_  windward  (Cp  )  is  read  from 
Graphs  D,  E,  and  F  depending  on  the  value  of  side  ratio.  The  pressure 
coefficient  Cpj  leeward  (Cp^)  is  taken  to  be  zero,  Cpj  =  0. 

3.  Calculate  the  local  wind  speed  at  mid-floor  height. 
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Table  B-2.  Summary  of  Wind  Pressure  Difference  Characteristics  at 
Mid-Height  of  Walls  on  Six  Types  of  Isolated  Models  at 
Wind  Incidences,  Normal,  30  Degrees  and  45  Degrees  to 
Long  Walls 


Model  Features 
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Single  story  on-grade 


Worst 


0 

30 

45 


0.49 

0.41 

0.34 


0.56  C 
0.67  WE 
0.71  WE 


0.42  BE 
0.10  LE 
-0.02  LE 


Average 


0.41 


Single  story  on-grade 
with  extended  eaves  and 
end  walls 

— 

2 

2nd 

Best 

0 

30 

45 

n 

0.56  C 
0.71  WE 
0.63  C 

0.46  BE 
0.58  LE 
0.52  C 

Average 

0.58 

Single  story  elevated 

3 

3rd 

0 

0.65 

0.71  C 

0.60  BE 

above  ground 

Best 

30 

0.58 

0.79  WE 

0.33  LE 

45 

0.43 

0.81  WE 

0.00  LE 

Average 

0.55 

Single  story  elevated 

4 

Best 

0 

0.85  BE 

0.77  C 

above  ground  with  extended 

30 

■in 

0.75  LE 

0.50  WE 

verandas  and  end  walls 

45 

B&B 

0.81  LE 

0.73  BE 

Average 

0.74 

Two  story 

5 

5  th 

0 

0.61 

0.69  C 

0.52  BE 

Best 

30 

0.33 

0.37  C 

0.24  WE 

45 

0.31 

0.42  C 

0.12  LE 

Average 

0.42 

Single  story  on-grade  with 

6 

4th 

0 

0.51 

0.58  C 

0.44  BE 

flat  roof 

Best 

30 

0.42 

0.69  WE 

0.08  LE 

45 

0.35 

0.71  WE 

0.00  LE 

Average 

0.43 

C  =  center  one-third  of  long  walls;  BE  =  both  ends  of  long  walls;  WE  =  windward  ends 
of  long  walls;  LE  =  leeward  ends  of  long  walls 
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W  =  smaller  building  width 
L  =  larger  building  width 


calculate  side  ratio:  y 


_ width  of  building 

length  (larger  side)  of  building 


Figure  B-2a.  Coordinate  system  of  building. 
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Figure  B-2b.  Wind  directions  on  buildings. 


Graphs:  Pressure  coefficients  for  buildings  three  stories  height  or  less  with  a  single 
value  reference  speed,  Ur,  at  building  roof  height,  ACp  =  Cp^  (Cp^  =  0). 
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A.  Averaged  pressure  coefficients  based  on  a  reference  velocity  measured 
at  the  roof,  Cp^,  for  a  side  ratio  of  1.0. 
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B.  Averaged  pressure  coefficients  based  on  a  reference  velocity  measured 
at  the  roof,  Cpp,  for  a  side  ratio  of  0.5. 
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Gamma  0.25 
O  Side  1 
A  Side  2 
□  Side  3 
O  Side  4 
O  Roof 


C  Averaged  pressure  coefficients  based  on  a  reference  velocity  measured 
at  the  roof,  Cpj>,  for  a  side  ratio  of  0.25. 


Graphs:  Pressure  coefficients  for  tall  buildings  for  multiple  values  of  reference 


velocity,  Uz,  ACp  =  Cpj. 


D.  Averaged  local  pressure  coefficients,  Cp^,  for  a  side  ratio  of  1.0. 
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E.  Averaged  local  pressure  coefficients,  Cp^,  for  a  side  ratio  of  0.5. 
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F.  Averaged  local  pressure  coefficients,  Cp^,  for  a  side  ratio  of  0.25. 


Appendix  C 

DISCHARGE  COEFFICIENTS  AND  AIR  LEAKAGE  DATA 
FOR  BUILDING  COMPONENTS 
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Table  C~2.  Summary  of  Air  Leakage  Data  for  Building  Components  Taken  From  Various  Sources 
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DISTRIBUTION  LIST 


ARMY  Fal  Engr.  Lelterkenny  Army  Depot,  Chamhersburg.  PA 
AF  HO  Traffic  Mgmnt  Cargo  Br  (LTC  S.  Rohrbough)  Washington. DC 

AFB  (AFITLDE).  Wright  Patterson  OH.  (RDVA)  AFESC  R&D  Tyndall.  FI..  X2ABC.  DEMC.  Williams  AZ: 
ABG/DEE  (F.  Nethers).  C.oodfellow  AFB  IX.  AF  Tech  Office  (Mgt  &  Ops).  Tyndall.  FL;  AFBSC  1ST. 
Library.  Tyndall.  FL:  CESCH.  Wrighl-Patterson:  Elec  Sys  Div.  Code  OCRS  (LTC  M.  Traister)  Hanseom. 
MA;  HO  MACDEEE.  Scott.  II:  SAMSO  MNND.  Norton  AFB  (  A:  Samso  Dee  (Sauer)  Vandenburg.  CA: 
Stinfo  Library.  Offutt  NE 
AFESC  DEB.  Tyndall.  FL 
NATL  ACADEMY  OF  ENG  Alexandria.  VA 

ARMY  ARRADCOM.  Dover.  NJ:  BMDSC-RE  (II.  McClellan)  Huntsville  AL:  C  ontracts  -  Faes  Engr 
Directorate.  Fort  Ord.  CA;  DAEN-CWE-M.  Washington  DC:  DAEN-MPE-D  Washington  DC: 

DAEN-MPU.  Washington  DC:  ERADCOM  Lech  Supp  Dir  (DELSD-I.)  Ft.  Monmouth.  NJ:  Natick  R&D 
Command  (Kwoh  Hu)  Natick  MA;  R&D  Command  DRDNA-l’STI  (J  Siegel)  Natick.  MA:  lech.  Ref. 

Div..  Fort  Huachuca  AZ 
ARMY  -  CERL  Library.  Champaign  IL 

ARMY  CORPS  OF  ENC.INEERSMRD-Eng  Div  .  Omaha  NE:  Seattle  D.st  Library .  Seattle  WA 

ARMY  CRREL  G.  Phetteplace  Hanover.  NH 

ARMY  DARCOM  Code  DRCMM-CS  Alexandria  VA 

ARMY  ENG  WATERWAYS  EXP  STA  Library.  Vicksburg  MS 

ARMY  ENGR  DIST.  Library.  Portland  OR 

ARMY  ENVIRON  HYGIENE  AGCY  HSE-EW  Water  Oual  Eng  Div  Aberdeen  Prov  Cirnd  MD 

ARMY  LOGISTICS  COMMAND  Code  ALCATCL-MS  (Morrissetl)  Fort  lee.  VA 

ARMY  MATERIALS  &  MECHANICS  RESEARCH  CENTER  Dr  l.cnoe.  W.i  ,-rtown  MA 

ARMY  MISSILE  R&D  CMD  SCI  Info  Cen  (DOC)  Redstone  Arsenal.  AL 

ARMY  MTMC  Trans  Engr  Agency  MTT-CE.  Newport  News.  VA 

ARMY  TRANSPORTATION  SCHOOL  Code  ATSPO  (  D  IE  Fort  Eustis.  VA 

ASO  PWD  (ENS  M  W  Davis).  Phildadclphia.  PA 

BUMED  Security  Offr.  Washington  DC 

BUREAU  OF  RECLAMATION  Code  1512  <C.  Selander)  Denver  CO 
CINCLANT  CIV  ENGR  SUPP  PLANS  OFFR  NORFOLK.  VA 
CINCPAC  Fac  Engrng  Div  (J44)  Makalapa.  HI 
CNAVRES  Code  13  (Dir.  Facilities)  New  Orleans.  LA 

C'NM  Code  MAT-04.  Washington.  DC:  Code  MAT-OKE.  Washington.  DC:  NMAT  -  044.  Washington  DC 
('NO  Code  NOP-%4.  Washington  DC:  Code  OP  <JS7  Washington  DC:  Code  OP-413  Wash.  DC;  Code  OPNAV 
09B24  (H):  OPW7J.  Washington.  DC 
COMCBPAC  Operations  Off.  Makalapa  HI 

COMFLEACT.  OKINAWA  PWD  -  Engr  Div.  Sasebo,  Japan:  PWO.  Kadena.  Okinawa:  PWO.  Sasebo.  Japan 

COMNAVBEACHPHIBREFTRAGRU  ONE  San  Diego  CA 

COMNAVMARIANAS  Code  N4.  Guam 

COMOCEANSYSLANT  PW-FAC  MGMNT  Off  Norfolk.  VA 

COMOCEANSYSPAC  SCE.  Pearl  Harbor  HI 

COMSUBDEVGRUONE  Operations  Offr.  San  Diego.  CA 

COMUSNAVCENT  N42.  Pearl  Harbor.  HI 

DEFENSE  INTELLIGENCE  AGENCY  DB-4C1  Washington  DC 

DEFFUELSUPPCEN  DFSC-OWE  (Term  Engrng)  Alexandria.  VA;  DFSC-OWE.  Alexandria  VA 
DLSIE  Army  Logistics  Mgt  Center.  Fort  Lee.  VA 

DOE  Div  Ocean  Energy  Sys  Cons/Solar  Energy  Wash  DC;  INEL  Tech.  Lib.  (Reports  Section).  Idaho  Falls.  ID 
DTIC  Defense  Technical  Info  Ctrl  Alexandria.  VA 
DTNSRDC  Code  4111  (R.  C.ierich).  Bethesda  MD 
DTNSRDC  Code  522  (Library).  Annapolis  MD 

ENVIRONMENTAL  PROTECTION  AGENCY  Reg.  Ill  Library.  Philadelphia  PA:  Reg.  VIII.  XM-ASL. 

Denver  CO 

FLTCOMBATTRACENLANT  PWO.  Virginia  Bch  VA 
FMFLANT  CEC  Offr.  Norfolk  VA 
GIDEP  OIC.  Corona.  CA 

GSA  Assist  Comm  Des  &  Cnst  (FAIA)  D  R  Dibner  Washington.  DC  :  Off  of  Des  &  Const-PCDP  (D  Eakin) 
Washington.  DC 

HC  &  RS  Tech  Pres.  Service.  Medcn.  Washington.  DC 
KWAJALEIN  MISRAN  BMDSC'-RKL-C 

LIBRARY  OF  CONGRESS  Washington.  DC  (Sciences  &  Tech  Div) 

MARCORPS  Code  LME-I  (Riggs)  Washington.  DC 


MARINE  CORPS  BASE  Code  406.  Camp  Lcjcunc.  NC':  Mainl  Off  Camp  Pendleton,  CA:  PWD  -  Maim. 
Control  Di\.  Camp  Butler.  Kawasaki.  Japan:  PWO  Camp  l.ejeune  NC:  PWO.  Camp  Pendleton  CA;  PWO. 
Camp  S.  D  Butler.  Kawasaki  Japan 
MARINE  CORPS  HOS  Code  LFF-2.  Washington  DC 

MCAS  Facil.  Engr.  Div,  Cherry  Point  NC:  CO.  Kaneohe  Bay  111:  Code  S4.  Ouantieo  VA:  Fats  Main!  Dept  - 
Operations  Div.  Cherry  Point:  PWD  -  Utilities  Div.  Iwakuni.  Japan:  PWO.  iwukuni.  Japan:  PWO.  Yuma 
AZ 

MC’DEC  M&L  Div  (LTC  R.  Kcft)  Ouantieo.  VA:  M&L  Div  Ouantieo  VA.  NSAP  REP.  Ouantieo  VA 
MCLB  B52JI.  Barstow  CA;  Maintenanee  Officer.  Barstow.  CA;  PWO.  Barstow  CA 
MCRD  SCE.  San  Diego  CA 

NAF  PWD  -  Engr  Div.  Atsugi.  Japan:  PWO.  Atsugi  Japan 
NALF  OINC.  San  Diego.  CA 

NARF  Code  100.  Cherry  Point.  NC.  Code  612.  Jax.  FI.:  Code  640.  Pensacola  FL:  SCE  Norfolk.  VA 
NAS  CO.  Guantanamo  Bay  Cuba:  Code  114.  Alameda  CA:  Code  183  (Fae.  Plan  BR  MGR):  Code  187. 
Jacksonville  FL:  Code  187(H).  Brunswick  ME:  Code  I8U  (ENS  P.J.  Hickey).  Corpus  Christi  TX;  Code  8E. 
Patuxent  Riv..  MD:  Dir  of  Engrng.  PWD.  Corpus  Christi.  TX:  Dir.  Util.  Div  .  Bermuda:  Grover.  PWD. 
Patuxent  River.  MD:  Lakehurst.  NJ;  Lead.  Chief.  Petty  Offr.  PW  Self  Help  Div.  Beevillc  TX.  PW  (J. 
Maguire).  Corpus  Christi  TX:  PWD  -  Engr  Div  Dir.  Millington.  TN:  PWD  -  Engr  Div.  Gtmo.  Cuba:  PWD  - 
Engr  Div.  Oak  Harbor.  WA:  PWD  Maint.  Cont.  Dir..  Fallon  NV;  PWD  Maim.  Div..  New  Orleans.  Belle 
Chasse  LA:  PWD.  Code  1 82 1 H  (Pfankuch)  Miramar.  SD  CA:  PWD.  Maintenance  Control  Dir..  Bermuda: 
PWD.  Willow  Grove  PA:  PWO  Belle  Chasse.  LA:  PWO  Chase  Field  Beevillc.  TX:  PWO  Key  West  FL: 
PWO  Lakehurst.  NJ:  PWO  Sigonella  Sicily:  PWO  Whiting  Fid.  Milton  FL:  PWO.  Dallas  TX:  PWO. 

Glenview  IL:  PWO,  Millington  TN:  PWO.  Miramar.  San  Diego  CA:  SCE  Norfolk.  VA:  SCE.  Barbers  Point 
HI:  SCE.  Cubi  Point.  R.P 
NASDC-WDC  T.  Fry.  Manassas  VA 

NATL  BUREAU  OF  STANDARDS  R  Chung  Washington.  DC 
NATL  RESEARCH  COUNCIL  Naval  Studies  Board.  Washington  DC 
NAVACT  PWO.  London  UK 
NAVAEROSPREGMEDCEN  SCE.  Pensacola  FL 

NAVAIRDEVCEN  Chmielewski.  Warminster.  PA:  PWD.  Engr  Div  Mgr.  Warminster.  PA 
NAVAIRPROPTESTCEN  CO.  Trenton.  NJ 
NAVAIRSYSCOM  Code  417I2A  (Michols)  Washington.  DC 
NAVAIRTESTCEN  PATUXENT  RIVER  PWD  (F.  McGrath).  Patuxent  Riv  .MD 
NAVCHAPGRU  Operations  Officer.  Code  3(1  Williamsburg.  VA 

NAVCOASTSYSCEN  CO.  Panama  City  FL;  Code  715  (J  Quirk)  Panama  City.  FL;  Code  772  (C  B  Kocsy) 
Panama  Citv  FL:  Librarv  Panama  City.  FL:  PWO  Panama  City.  FL 
NAVC'OMMAREAMSTRSTA  PWO.  Norfolk  VA;  SCE  Unit  I  Naples  Italy:  SCE.  Wahiawa  HI 
NAVCOMMSTA  Code  401  Nea  Makri.  Greece:  PWD  -  Maint  Control  Div.  Diego  Garcia  Is.:  PWO.  Exmouth. 
Australia:  SCE.  Balboa.  CZ 

NAVCONSTRACEN  Code  00UI5,  Port  Hueneme  CA:  Curriculum/lnstr.  Stds  Offr.  Gulfport  MS 
NAVEDTRAPRODEVCEN  Technical  Library.  Pensacola.  FL 
NAVEDUTRACEN  Engr  Dept  (Code  42)  Newport.  Rl 
NAVEODTECHCEN  Code  60.5.  Indian  Head  MD 

NAVFAC  PWO.  Brawdy  Wales  UK:  PWO.  Centerville  Bch.  Ferndale  CA:  PWO.  Point  Sur.  Big  Sur  CA 
NAVFACENGCOM  Alexandria.  VA;  Code  03  Alexandria.  VA:  Code  03T  (Essoglou)  Alexandria.  VA:  Code 
043  Alexandria.  VA;  Code  044  Alexandria.  VA:  Code  044B)  Alexandria.  VA:  Code  0451  (P  W  Brewer) 
Alexandria.  Va;  Code  0451,  Alexandria.  VA;  Code  (MAI  Alexandria.  VA;  Code  (MB3  Alexandria.  VA: 

Code  05IA  Alexandria,  VA:  Code  06.  Alexandria  VA;  Code  OVM54.  Tech  Lib.  Alexandria.  VA:  Code  1(H) 
Alexandria.  VA;  Code  1113.  Alexandria.  VA:  Code  1 1  IB  Alexandria.  VA:  code  081  Alexandria.  VA 
NAVFACENGCOM  -  CHES  DIV.  Code  403  Washington  DC:  Code  405  Wash.  DC;  FPO-I  Washington.  DC: 
Library.  Washington.  D  C. 

NAVFACENGCOM  -  LANT  DIV.  Code  III.  Norfolk.  VA;  Code  403.  Norfolk.  VA;  Code  405  Civil  Engr  BR 
Norfolk  VA:  Eur.  BR  Deputy  Dir.  Naples  Italy;  Library.  Norfolk.  VA;  RDT&ELO  I02A.  Norfolk.  VA 
NAVFACENGCOM  -  NORTH  DIV.  Code  (M  Philadelphia.  PA:  Code  OOP  Philadelphia  PA:  Code  1028. 
RDT&ELO.  Philadelphia  PA;  Code  III  Philadelphia.  PA.  Library.  Philadelphia.  PA:  ROICC.  Contracts. 
Crane  IN 

NAVFACENGCOM  -  PAC  DIV.  (Kyi)  Code  101.  Pearl  Harbor.  HI:  CODE  OOP  PEARL  HARBOR  HI:  Code 
2011  Pearl  Harbor.  HI:  Code  402.  RDT&E.  Pearl  Harbor  HI;  Commander.  Pearl  Harbor.  HI:  Library. 

Pearl  Harbor.  HI 

NAVFACENGCOM  -  SOUTH  DIV.  Code  403.  Gaddy.  Charleston.  SC:  Code  00.  RDT&ELO.  Charleston  SC: 
Library.  Charleston,  SC 

NAVFACENGCOM  -  WEST  DIV.  AROICC.  Contracts.  Twenty-nine  Palms  CA.  Code  IMB  San  Bruno.  CA: 

Library.  San  Bruno.  CA:  OOP/21)  San  Bruno.  CA:  RDT&ELO  Code  2011  San  Bruno.  CA 
NAVFACENGCOM  CONTRACTS  AROICC.  NAVSTA  Brooklyn.  NY:  AROICC.  Ouantieo.  VA;  C  ontracts. 
AROICC.  Lemoorc  CA:  Dir.  Eng.  Div..  Exmouth.  Australia:  Eng  Div  dir.  Southwest  Pac.  Manila.  PI: 

OICC.  Southwest  Pac.  Manila.  PI:  OICC-ROIC'C.  NAS  Oceana.  Virginia  Beach.  VA:  OICC  ROICC. 


Balboa  Panama  Canal:  ROICC  AF  Guam:  ROK'C  Code  495  Portsmouth  VA:  ROICC  Key  West  FI.: 

ROICC  MCAS  El  Toro:  ROICC.  Keflavik.  Iceland:  ROICC.  NAS.  Corpus  Christi.  IX:  ROICC  .  Pacific.  San 
Bruno  CA:  ROICC.  Yap:  ROICC-OICC-SPA.  Norfolk.  VA 
NAVHOSP  PWD  -  Engr  Div.  Beaufort.  SC 
NAVMAG  PWD  •  Engr  Div.  Guam;  SCE.  Subie  Bay.  R.P. 

NAVOCEANSYSCEN  Code  4473B  (Tech  Lib)  San  Diego.  CA.  Code  523  (Hurley).  San  Diego.  CA.  Code  67tHt. 

San  Diego.  CA:  Code  SI  I  San  Diego.  CA 
NAVORDMISTESTFAC  PWD  -  Engr  Dir.  White  Sands.  NM 
NAVORDSTA  PWD  -  Dir.  Engr  Div.  Indian  Head.  MD;  PWO.  Louisville  KY 
NAVPETOFF  Code  30.  Alexandria  VA 
NAVPETRES  Director.  Washington  DC 

NAVPHIBASE  CO.  AC'B  2  Norfolk.  VA:  Code  S3T.  Norfolk  VA.  SCE  Coronado.  SD.CA 
NAVREGMEDCEN  PWD  -  Engr  Div.  Camp  Lejeune.  NC:  PWO.  Camp  Lejeunc.  NC 
NAVREGMEDCEN  PWO.  Okinawa.  Japan 

NAVREGMEDCEN  SCE;  SCE  San  Diego.  CA:  SCE.  Camp  Pendleton  CA:  SCE.  Guam:  SCE.  Newport.  Rl: 
SCE.  Oakland  CA 

NAVREGMEDCEN  SCE.  Yokosuka.  Japan 
NAVSCOLCECOFF  C'35  Port  Hueneme.  CA 
NAVSCSOL  PWO.  Athens  GA 

NAVSEASYSCOM  Code  0325,  Program  Mgr.  Washington.  DC:  Code  PMS  305  A  3.  Washington.  DC:  SEA 
04E  (L  Kess)  Washington.  DC 

NAVSECGRUACT  PWO.  Adak  AK:  PWO.  Edzell  Scotland:  PWO.  Puerto  Rico:  PWO.  lorri  Sta.  Okinawa 
NAVSECSTA  PWD  -  Engr  Div.  Wash..  DC 

NAVSHIPYD  Code  202.4.  Long  Beach  CA;  Code  202.5  (Library)  Puget  Sound.  Bremerton  WA:  Code  3X0. 
Portsmouth.  VA:  Code  382.3.  Pearl  Harbor.  HI:  Code  400.  Puget  Sound:  Code  410.  Mare  Is..  Vail-.' jo  CA. 
Code  440  Portsmouth  NH;  Code  440,  Norfolk;  Code  440.  Puget  Sound.  Bremerton  WA:  Code  453  (I'til. 
Supr).  Vallejo  CA:  Library.  Portsmouth  NH.  PW  Dept.  Long  Beach.  CA;  PWD  (Code  420)  Dir  Portsmouth. 
VA.  PWD  (Code  450-HD)  Portsmouth.  VA:  PWD  (Code  453-HD)  SHPO  03.  Portsmouth.  VA:  PWO. 
Bremerton.  WA:  PWO.  Mare  Is.:  PWO.  Puget  Sound:  SCE.  Pearl  Harbor  HI;  Tech  Library.  Vallejo.  CA 
NAVSTA  Adak.  AK:  CO.  Brooklyn  NY:  Code  4.  12  Marine  Corps  Dist.  Treasure  Is..  San  Francisco  CA;  Dir 
Engr  Div.  PWD.  Mavport  FL:  Dir  Mech  Engr  37WC93  Norfolk.  VA:  Engr.  Dir..  Rota  Spain:  Long  Beach. 
CA:  Maint.  Corn.  Div..  Guantanamo  Bay  Cuba:  Maint.  Div.  Dir  Code  531.  Rodman  Panama  Canal:  PWD  • 
Engr  Dept.  Adak.  AK:  PWD  -  Engr  Div.  Midway  Is.:  PWO.  Keflavik  Iceland:  PWO.  Mavport  FL:  SCE. 
Guam:  SCE.  Pearl  Harbor  HI:  SCE.  San  Diego  CA:  SCE.  Subic  Bay.  R.P.:  Utilities  Engr  Off.  Rota  Spain 
NAVSUPPACT  CO.  Naples.  Italy:  PWO  Naples  Italy 
NAVSUPPFAC  PWD  -  Maint  Control  Div.  Thurmont.  MD 
NAVSURFWPNCEN  PWO.  White  Oak.  Silver  Spring.  MD 
NAVTECHTRACEN  SCE.  Pensacola  FL 
NAVTELCOMMCOM  Code  53.  Washington.  DC 

NAVWPNCEN  Code  2636  China  Lake:  Code  3803  China  Lake.  CA:  PWO  (Code  266)  China  Lake.  CA:  ROUT 
(Code  702).  China  Lake  CA 

NAVWPNSTA  (Clebak)  Colts  Neck.  NJ;  Code  002.  Concord  CA:  Code  002 A.  Seal  Beach.  CA 
NAVWPNSTA  PW  Office  Yorktown.  VA 

NAVWPNSTA  PWD  -  Maint.  Control  Div..  Concord.  CA;  PWD  -  Supr  Gen  Engr.  Seal  Beach.  CA.  PWO. 

Charleston.  SC:  PWO.  Seal  Beach  CA 
NAVWPNSUPPCEN  Code  09  Crane  IN 
NCTC  Const.  Elec.  School.  Port  Hueneme.  CA 

NCBC  Code  10  Davisville.  Rl:  Code  15.  Port  Hueneme  CA;  Code  155.  Port  Hueneme  CA.  Code  1552  (Bra/ele) 
Port  Hueneme.  CA:  Code  I55B  (Nishimura)  Port  Hueneme.  CA:  Code  156.  Port  Hueneme.  CA:  Code  I56F 
(Volpe)  Port  Hueneme.  CA:  Code  251 1 1  Port  Hueneme.  CA:  Code  430  (PW  Engrng)  Gulfport.  MS:  Code 
470.2.  Gulfport.  MS.  NEESA  Code  252  (P  Winters)  Port  Hueneme.  CA:  PWO  (Code  80)  Port  Hueneme. 

CA:  PWO.  Davisville  Rl:  PWO.  Gulfport.  MS 
NCR  20.  Code  R70;  20.  Commander 

NMCB  3.  SWC  D.  Wellington;  74.  CO:  FIVE.  Operations  Dept;  Forty,  CO:  THREE.  Operations  Off 

NOAA  (Dr.  T.  Me  Guinness)  Rockville.  MD:  Library  Rockville.  MD 

NRL  Code  5800  Washington.  DC 

NROTC  J.W.  Stephenson.  UC.  Berkeley.  CA 

NSC  Code  54.1  Norfolk.  VA 

NSD  SCE.  Subic  Bay.  R.P. 

NSWSES  Code  0150  Port  Hueneme.  CA 
NTC  OICC.  CBU-401.  Great  Lakes  IL 

NUSC  Code  131  New  London.  CT:  Code  5202  (S.  Schady)  New  London.  CT:  Code  EAI23  (R.S  Munn).  New 
London  CT:  Code  SB  331  (Brown).  Newport  Rl 

OFFICE  SECRETARY  OF  DEFENSE  OASD  (MRA&L)  Dir.  of  Energy.  Pentagon.  Washington.  DC 
ONR  Code  221.  Arlington  VA:  Code  700F  Arlington  VA 
PACMISRANFAC'  HI  Area  Bkg  Sands.  PWO  Kekaha.  Kauai.  HI 


PHIBCB  I  P&E.  San  Diego.  (  .»:  I.  CO  San  Diego.  (  A 

PWC  ACE  Office  Norfolk.  VA  CO  Norfolk.  VA:  CO.  (Code  10).  Oakland.  CA:  CO.  Great  l  akes  II.:  CO. 
Pearl  Harbor  HI:  Code  10.  Great  lakes.  II.;  Code  105  Oakland.  CA:  Code  III).  Great  Lakes.  II.:  Code  110. 
Oakland.  CA:  Code  120.  Oakland  CA:  Code  I2K.  Guam:  Code  154  (Library).  Great  Lakes.  II.:  Code  200. 
Great  Lakes  IL:  Code  400.  Great  Lakes.  II.;  Code  400.  Oakland.  CA:  Code  400.  Pearl  Harbor.  HI;  Code 
400.  San  Diego.  CA:  Code  420.  Great  Lakes.  II. ;  Code  420.  Oakland.  CA:  Code  424.  Norfolk.  VA:  Code 
500  Norfolk.  VA:  Code  505A  Oakland.  CA:  Code  (410.  Great  Lakes.  IL:  Code  bill.  San  Diego  Ca:  Code 
700.  Great  Lak -s.  IL:  Library.  Code  I20C.  San  Diego.  CA;  Library.  Guam:  Library.  Norfolk.  VA:  l.ihrarv. 
Oakland.  C'A:  Library.  Pearl  Harbor.  HI:  Library.  Pensacola.  EL;  Library.  Subic  Bay.  R  P  :  Library. 
Yokosuka  JA:  Util  Dept  (R  Paseua)  Pearl  Harbor.  HI:  Utilities  Officer.  Guam 
SPCC  PWO  (Code  120)  Mechaniesburg  PA 
SUPANX  PWO.  Williamsburg  VA 

TVA  Smelser.  Knoxville.  Tenn.;  Solar  Group.  Arnold.  Knoxville.  IN 
U  S.  MERCHANT  MARINE  ACADEMY  Kings  Point.  NY  (Reprint  Custodian) 

US  NAVAL  FORCES  Korea  (ENJ-P&O) 

USAF  REGIONAL  HOSPITAL  Fairchild  AFB.  WA 
IISCG  (Smith).  Washington.  DC:  G-MMT-482  (J  Spencer) 

USDA  Forest  Service  Reg  3  (R.  Brown)  Albuquerque.  NM 

USNA  Ch.  Mcch.  Engr.  Dept  Annapolis  MD:  ENCiRNG  Div.  PWD.  Annapolis  MD:  Energy- Environ  Study 
Grp.  Annapolis.  MD.  Environ.  Prot.  R&D  Prog.  (J.  Williams).  Annapolis  MD:  Mech.  Engr.  Dept.  (C. 

Wu).  Annapolis  MD:  USNA/Sys  Eng  Dept.  Annapolis.  MD 
USS  FULTON  WPNS  Rep.  Offr  (W-3)  New  York.  NY 
ARIZONA  Kroelinger  Tempe.  AZ:  State  Energy  Programs  Off..  Phoenix  AZ 
AUBURN  IINIV.  Bldg  Sci  Dept.  Lechner.  Auburn.  AL 
BERKELEY  PW  Engr  Div.  Harrison.  Berkeley.  CA 

BONNEVILLE  POWER  ADMIN  Portland  OR  (Energy  Consrv.  Off..  D  Davey) 

BROOKHAVEN  NATL  LAB  M.  Steinberg.  Upton  NY 

CALIFORNIA  STATE  UNIVERSITY  LONG  BEACH.  CA  (CHELAPATI) 

CONNECTICUT  Office  of  Policy  &  Mgt.  Energy.  Div.  Hartford.  CT 
CORNELL  UNIVERSITY  Ithaca  NY  (Serials  Dept.  Engr  Lib.) 

DAMES  &  MOORE  LIBRARY  LOS  ANGELES.  CA 

DRURY  COLLEGE  Physics  Dept.  Springfield.  MO 

FLORIDA  ATLANTIC  UNIVERSITY  Boca  Raton.  FL  (McAllister) 

FOREST  INST  FOR  OCEAN  &  MOUNTAIN  Carson  City  NV  (Studies  -  Library) 

GEORGIA  INSTITUTE  OF  TECHNOLOGY  (LT  R.  Johnson)  Atlanta.  GA:  Col.  Arch.  Benton.  Atlanta.  GA 

HARVARD  UNIV.  Dept,  of  Architecture.  Dr.  Kim.  Cambridge.  MA 

HAWAII  STATE  DEPT  OF  PLAN  &  ECON  DEV.  Honolulu  HI  (Tech  Info  Ctrl 

IOWA  ST  ATE  UNIVERSITY  Dept  Arch.  McKrown.  Ames.  IA 

WOODS  HOLE  OCEANOGRAPHIC  INST  Woods  Hole  MA  (Winget) 

KEENE  STATE  COLLEGE  Keene  NH  (Cunningham) 

LEHIGH  UNIVERSITY  Bethlehem  PA  (Linderman  Lib.  No. 30.  Flecksteiner) 

LOUISIANA  DIV  NATURAL  RESOURCES  &  ENERGY  Div  Of  R&D.  Baton  Rouge.  LA 
MAINE  OFFICE  OF  ENERGY  RESOURCES  Augusta.  ME 
MISSOURI  ENERGY  AGENCY  Jefferson  City  MO 

MIT  Cambridge  MA:  Cambridge  MA  (Rm  10-500.  Tech.  Reports.  Engl.  Lib.):  Cambridge.  MA  (Harlcman) 

MONTANA  ENERGY  OFFICE  Anderson.  Helena.  MT 

NATL  ACADEMY  OF  ENG.  ALEXANDRIA.  VA  (SEARLE.  JR.) 

NATURAL  ENERGY  LAB  Library.  Honolulu.  HI 

NEW  HAMPSHIRE  Concord  NH  (Governor's  Council  on  Energy) 

NEW  MEXICO  SOLAR  ENERGY  INST.  Dr.  Zwibcl  Las  Cruces  NM 
NY  CITY  COMMUNITY  COLLEGE  BROOKLYN.  NY  (LIBRARY) 

NYS  ENERGY  OFFICE  Library.  Albany  NY 

OAK  RIDGE  NATL  LAB  T.  Lundy.  Oak  Ridge.  TN 

PURDUE  UNIVERSITY  Lafayette.  IN  <CE  Engr  Lib) 

SCRIPPS  INSTITUTE  OF  OCEANOGRAPHY  LA  JOLLA.  CA  (ADAMS) 

SEATTLE  U  Prof  Schwaegler  Seattle  WA 

SRI  INTL  Phillips.  Chcm  Engr  Lab.  Menlo  Park.  CA 

STATE  UNIV.  OF  NEW  YORK  Fort  Schuyler.  NY  (Longobardi) 

STATE  UNIV  OF  NY  AT  BUFFALO  School  of  Medicine.  Buffalo.  NY 
TEXAS  A&M  UNIVERSITY  W.B.  Ledbetter  College  Station.  TX 

UNIVERSITY  OF  CALIFORNIA  Energy  Engineer.  Davis  CA:  LIVERMORE.  CA  (LAWRENCE 
LIVERMORE  LAB.  TOKARZ);  UCSF.  Physical  Plant.  San  Francisco.  CA 
UNIVERSITY  OF  DELAWARE  Newark.  DE  (Dept  of  Civil  Engineering.  Chcsson! 

UNIVERSITY  OF  FLORIDA  Dept  Arch..  Morgan.  Gainesville.  FL 
UNIVERSITY  OF  HAWAII  HONOLULU.  HI  (SCIENCE  AND  TECH  DIV.) 

UNIVERSITY  OF  ILLINOIS  (Hall)  Urbana.  IL;  Metz  Ref  Rm.  Urbana  IL:  URBANA.  IL  (LIBRARY) 


UNIVERSITY  OF  MASSACHUSETTS  (Heronemus).  ME  Dept  Amherst.  MA 
UNIVERSITY  OF  NEBRASKA-LINCOLN  Lincoln.  NE  (Ross  Icc  Shell  Proj.) 

UNIVERSITY  OF  NEW  HAMPSHIRE  Elec  Engr  Depot.  I)r  Murdoch.  Durham.  N  II 
UNIVERSITY  OF  TEXAS  Inst.  Marine  Set  (Library!.  Port  Arkansas  TX 
UNIVERSITY  OF  TEXAS  AT  AUSTIN  AUSTIN.  IX  (THOMPSON) 

UNIVERSITY  OF  WASHINGTON  Seattle  WA  (E.  Linger) 

UNIVERSITY  OF  WISCONSIN  Milwaukee  Wl  (Ctr  of  Great  Lakes  Studies) 

ANALYTICAL  TECH  Lawrence.  Port  Hucncme.  CA 
ARVID  GRANT  OLYMPIA.  WA 
ATLANTIC  RICHFIELD  CO.  DALLAS.  TX  (SMI  I  II  I 
AWWA  RSCH  FOUNDATION  R  Heaton.  Denver  CO 
BECHTEL  CORP  SAN  FRANCISCO.  CA  (PHELPS) 

BROWN  &  ROOT  Houston  TX  (D.  Ward) 

CHEMED  CORP  Lake  Zurich  IL  (Dearborn  t  hem  Div  Lib  ) 

COLUMBIA  GULF  TRANSMISSION  CO.  HOUSTON.  TX  (ENG.  LIB.) 

DESIGN  SERVICES  Beck.  Ventura.  CA 

DIXIE  DIVING  CENTER  Decatur.  GA 

DURLACH.  O  NEAL.  JENKINS  &  ASSOC.  Columbia  SC 

FURGO  INC.  Library.  Houston.  TX 

CARD  INC.  Dr.  L.  Holmes.  Niles.  IL 

LINEFAST  CORP  (J.  DiMartino)  Holbrook.  NY 

LITHONIA  LIGHTING  Application  eng.  Dept.  (B.  Helton).  Conyers.  GA  302(17 
MCDONNEL  AIRCRAFT  CO.  (Favman)  Engrng  Dept..  St.  Louis.  MO 
MEDERMOTT  &  CO.  Diving  Division.  Harvey.  LA 

NEWPORT  NEWS  SHIPBLDG  &  DRYDOCK  CO.  Newport  News  VA  (  lech.  Lib  ) 

PACIFIC  MARINE  TECHNOLOGY  (M.  Wagner)  Duvall.  WA 
PG&E  Library.  San  Francisco.  CA 

PORTLAND  CEMENT  ASSOC  SKOKIE.  IL  (CORLEY:  Skokie  IL  (Rsch  &  Dev  Lab.  Lib.) 
RAYMOND  INTERNATIONAL  INC  E  Colic  Soil  Tech  Dept.  Pennsauken.  NJ 
ROHR  IND  INC  (D.  Buchanan)  Chula  Vista.  CA 

SANDIA  LABORATORIES  Albuquerque.  NM  (Vortman):  Library  Div..  Livermore  CA 
SCHUPACK  ASSOC  SO.  NORWALK.  CT  (SC'HUPACK) 

SHELL  DEVELOPMENT  CO.  Houston  TX  (C.  Sellars  Jr.) 

TEXTRON  INC  BUFFALO.  NY  (RESEARCH  CENTER  LIB.) 

TRW  SYSTEMS  REDONDO  BEACH.  CA  (DAI) 

UNITED  TECHNOLOGIES  Windsor  Locks  CT  (Hamilton  Std  Div..  Librarv) 

WARD.  WOLSTENHOLD  ARCHITECTS  Sacramento.  CA 

WESTINGHOUSE  ELECTRIC  CORP.  Annapolis  MD  (Oceanic  Div  Lib.  Brvan):  Library.  Pittsburgh  PA 
WM  CLAPP  LABS  -  BATTELLE  DUXBURY.  MA  (LIBRARY) 

BRAHTZ  La  Jolla.  CA 
BULLOCK  La  Canada 
FISHER  San  Diego.  Ca 
KETRON.  BOB  Ft  Worth.  TX 
KRUZIC.  T.P.  Silver  Spring.  MD 
CAPT  MURPHY  Sunnyvale.  CA 
R.F.  BESIER  Old  Saybrook  CT 

BROWN  &  CALDWELL  Saunders.  EM  ./Oakland.  CA 
T.W.  MERMEL  Washington  DC 
WALTZ  Livermore.  CA 


